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Hadrons are composite systems of quarks and gluons. A direct consequence is the possibility to have 
hadron-hadron collisions in which several distinct pairs of partons collide with each other: multiple 
interactions, a.k.a. multiple scatterings. At first glance, the divergence of the perturbative t-channel one- 
gluon-exchange graphs in the p± — > limit implies an infinity of interactions per event. However, the 
perturbative framework does not take into account screening from the fact that a hadron is in an overall 
colour singlet state. Therefore an effective cutoff p± m - m of the order of one to a few GeV is introduced, 
representing an inverse colour correlation distance inside the hadron. For realistic pj_ m i n values most 
inelastic events in high-energy hadronic collisions should then contain several perturbatively calculable 



O '. interactions, in addition to whatever nonperturbative phenomena may be present. 

Although most of this activity is not hard enough to play a significant role in the description of 
""■^ high-px jet physics, it can be responsible for a large fraction of the total multiplicity (and large fluctu- 

ations in it), for semi-hard (mini-)jets in the event, for the details of jet profiles and for the jet pedestal 
effect, leading to random as well as systematic shifts in the jet energy scale. Thus, a good understanding 
[ of multiple interactions would seem prerequisite to carrying out precision studies involving jets and/or 

the underlying event in hadronic collisions. 

In an earlier study it was argued that all the underlying event activity is triggered by the 
multiple interactions mechanism. However, while the origin of underlying events is thus assumed to be 
perturbative, many nonperturbative aspects still need to be considered and understood: 

( i) What is the detailed mechanism and functional form of the dampening of the perturbative cross section 
at small p±? (Certainly a smooth dampening is more realistic than a sharp p± m i n cutoff.) 

( ii) Which energy dependence would this mechanism have? 

(Hi) How is the internal structure of the proton reflected in an impact-parameter-dependent multiple 
interactions rate, as manifested e.g. in jet pedestal effects? 

(iv) How can the set of colliding partons from a hadron be described in terms of correlated multiparton 
distribution functions of flavours and longitudinal momenta? 

(v) How does a set of initial partons at some low perturbative cutoff scale, 'initiators', evolve into such a 
set of colliding partons? (Two colliding partons could well have a common initiator.) Is standard DGLAP 
evolution sufficient, or must BFKL/CCFM effects be taken into account? 

(vi) How would the set of initiators correlate with the flavour content of, and the longitudinal momentum 
sharing inside, the left-behind beam remnant? 

(vii) How are the initiator and remnant partons correlated by confinement effects ('primordial k± , )7 

(viii) How are all produced partons, both the interacting and the beam-remnant ones, correlated in colour? 
Is the large number-of-colours limit relevant, wherein partons can be hooked up into strings (with quarks 
as endpoints and gluons as intermediate kinks) representing a linear confinement force [2 |? 

(ix) How is the original baryon number of an incoming proton reflected in the colour topology? 
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(x) To what extent would a framework with independently fragmenting string systems, as defined from 
the colour topology, be modified by the space-time overlap of several strings? 

Needless to say, we should not expect to find a perfect solution to any of these issues, but only 
successively improved approximations. The framework in Q is very primitive in a number of respects. 
Nevertheless, it has turned out to be quite successful. Thus the Pythia Tune A of R.D. Field is 
capable of describing a host of jet and minimum-bias event data at the Tevatron. The model appears 
inadequate to fully describe correlations and fluctuations, however, and we would expect a poor perfor- 
mance for several topics not yet studied experimentally. 

In particular, only very simple beam remnant structures could technically be dealt with in [ 1 1. One 
recent development was the extension of the standard Lund string framework |2| to include a junction 
fragmentation description |4j that allows the hadronization of nontrivial colour topologies containing 
non-zero baryon number. In the context of multiple interactions, this improvement means that almost 
arbitrarily complicated baryon beam remnants may now be dealt with, hence many of the restrictions 
present in the old model are no longer necessary. 

Here, we report on the development of a new model for the flavour-, colour-, and momentum- 
correlated partonic structure involved in a hadron-hadron collision, i.e. partly addressing several of the 
points above. We first present the main work on flavour and momentum space correlations, and thereafter 
separately the very thorny issue of colour correlations, before concluding. A more complete description 
of the model, also including references to experimental data and other theoretical ideas, and with com- 
ments on the all the issues, may be found in 0. A toy model study of the first two points is found in [ 6 1. 
The Pythia manual [ 7 ] contains some complementary information. 

2. CORRELATED PARTON DENSITIES 

Consider a hadron undergoing multiple interactions in a collision. Such an object should be described 
by multi-parton densities, giving the joint probability of simultaneously finding n partons with flavours 
fi, . . . , f n , carrying momentum fractions x±, . . . ,x n inside the hadron, when probed by interactions at 
scales Qi , . . . , Q 2 n . However, we are nowhere near having sufficient experimental information to pin 
down such distributions. Therefore, and wishing to make maximal use of the information that we do 
have, namely the standard one-parton-inclusive partem densities, we propose the following strategy. 

As described in Q, the interactions may be generated in an ordered sequence of falling p±. For 
the hardest interaction, all smaller p± scales may be effectively integrated out of the (unknown) fully 
correlated distributions, leaving an object described by the standard one-parton distributions, by defini- 
tion. For the second and subsequent interactions, again all lower-p^ scales can be integrated out, but the 
correlations with the first cannot, and so on. Thus, we introduce modified parton densities, that correlate 
the i'th interaction and its shower evolution to what happened in the i — 1 previous ones. 

The first and most trivial observation is that each interaction i removes a momentum fraction Xi 
from the hadron remnant. Already in (Q this momentum loss was taken into account by assuming a 
simple scaling ansatz for the parton distributions, f(x) — > f(x/X)/X, where X = 1 — Yh=i x i * s the 
momentum remaining in the beam hadron after the n first interactions. Effectively, the PDF's are simply 
'squeezed' into the range x G [0, X]. 

Next, for a given baryon, the valence distribution of flavour / after n interactions, qf vn (x, Q 2 ), 
should integrate to the number Nf vn of valence quarks of flavour / remaining in the hadron remnant. 
This rule may be enforced by scaling the original distribution down, by the ratio of remaining to original 
valence quarks Nf vn /Nf v Q, in addition to the x scaling mentioned above. 

Also, when a sea quark is knocked out of a hadron, it must leave behind a corresponding antisea 
parton in the beam remnant. We call this a companion quark. In the perturbative approximation the sea 
quark q s and its companion q c come from a gluon branching g — > q s + q c (it is implicit that if q s is a 
quark, q c is its antiquark). Starting from this perturbative ansatz, and neglecting other interactions and 



any subsequent perturbative evolution of the q c , we obtain the q c distribution from the probability that a 
sea quark q s , carrying a momentum fraction x s , is produced by the branching of a gluon with momentum 
fraction y, so that the companion has a momentum fraction x = y — x s , 



q c (x; xs) oc J g{y) P g ^q s q c (z) 5(x s - zy) dz = -Pg-qsqc (^^) ' (1) 

with P g wq s q c the usual DGLAP gluon splitting kernel. A simple ansatz g(x) oc (1 — x) n jx is here used 
for the gluon. Normalizations are fixed so that a sea quark has exactly one companion. Qualitatively, 
xq c (x; x s ) is peaked around x « x s , by virtue of the symmetric P g ->q s q c splitting kernel. 

Without any further change, the reduction of the valence distributions and the introduction of com- 
panion distributions, in the manner described above, would result in a violation of the total momentum 
sum rule, that the x- weighted parton densities should integrate to X: by removing a valence quark from 
the parton distributions we also remove a total amount of momentum corresponding to (xf v ), the average 
momentum fraction carried by a valence quark of flavour /, 

/ v _ So xq fvn (x,Q 2 ) dx 

( X fvn) = % —. = X (z/vo) , (2) 

Jo Qfvn{x,Q 2 ) dx 

and by adding a companion distribution we add an analogously defined momentum fraction. 

To ensure that the momentum sum rule is still respected, we assume that the sea+gluon normal- 
izations fluctuate up when a valence distribution is reduced and down when a companion distribution is 
added, by a multiplicative factor 

1 -J2fN fvn (x fv0 } - E/,j(g/cjo) (3) 
l-Y,f N fvo(xfvo) 

The requirement of a physical x range is of course still maintained by 'squeezing' all distributions into 
the interval x € [0, X}. The full parton distributions after n interactions thus take the forms 



M 2 ) = ^ 



Nf vn ( x 2 \ fx 



' 1 X.c 



(4) 



g n {x) = —g {-^>Q J ' ( 5 ) 

where qf v0 (g/ s o) denotes the original valence (sea) distribution of flavour /, and the index j on the 
companion distributions qf c . counts different companion quarks of the same flavour /. 

After the perturbative interactions have taken each their fraction of longitudinal momentum, the 
remaining momentum is to be shared between the beam remnant partons. Here, valence quarks receive 
an x picked at random according to a small-Q 2 valence-like parton density, while sea quarks must be 
companions of one of the initiator quarks, and hence should have an x picked according to the q c [x\ x s ) 
distribution introduced above. In the rare case that no valence quarks remain and no sea quarks need 
be added for flavour conservation, the beam remnant is represented by a gluon, carrying all of the beam 
remnant longitudinal momentum. 

Further aspects of the model include the possible formation of composite objects in the beam 
remnants (e.g. diquarks) and the addition of non-zero primordial k± values to the parton shower initiators. 
Especially the latter introduces some complications, to obtain consistent kinematics. Details on these 
aspects are presented in [5 1. 
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Fig. 1: (a) The initial state of a baryon, with the valence quarks colour-connected via a central string junction J. (b) Example of 
a topology with initiators connected at random, (c) Alternative with the junction in the remnant. 

3. COLOUR CORRELATIONS 

The initial state of a baryon may be represented by three valence quarks, connected antisymmetrically 
in colour via a central junction, which acts as a switchyard for the colour flow and carries the net baryon 
number, Fig.[^. 

The colour-space evolution of this state into the initiator and remnant partons actually found in 
a given event is not predicted by perturbation theory, but is crucial in determining how the system 
hadronizes; in the Lund string model [2], two colour-connected final state partons together define a 
string piece, which hadronizes by successive non-perturbative breakups along the string. Thus, the 
colour flow of an event determines the topology of the hadronizing strings, and consequently where 
and how many hadrons will be produced. The question can essentially be reduced to one of choosing 
a fictitious sequence of gluon emissions off the initial valence topology, since sea quarks together with 
their companion partners are associated with parent gluons, by construction. 

The simplest solution is to assume that gluons are attached to the initial quark lines in a random 
order, see Fig.m>. If so, the junction would rarely be colour-connected directly to two valence quarks in 
the beam remnant, and the initial-state baryon number would be able to migrate to large p± and small 
xf values. While such a mechanism should be present, there are reasons to believe that a purely random 
attachment exaggerates the migration effects. Hence a free parameter is introduced to suppress gluon 
attachments onto colour lines that lie entirely within the remnant, so that topologies such as Fig. QJ; 
become more likely. 

This still does not determine the order in which gluons are attached to the colour line between 
a valence quark and the junction. We consider a few different possibilities: 1) random, 2) gluons are 
ordered according to the rapidity of the hard scattering subsystem they are associated with, and 3) gluons 
are ordered so as to give rise to the smallest possible total string lengths in the final state. The two latter 
possibilities correspond to a tendency of nature to minimize the total potential energy of the system, 
i.e. the string length. Empirically such a tendency among the strings formed by multiple interactions is 
supported e.g. by the observed rapid increase of (p±) with n c h ar ged- It appears, however, that a string 
minimization in the initial state is not enough, and that also the colours inside the initial-state cascades 
and hard interactions may be nontrivially correlated. These studies are still ongoing, and represent the 
major open issues in the new model. 



4. CONCLUSION 

A new model for the underlying event in hadron-hadron collisions [ 5 ] has been introduced. This model 
extends the multiple interactions mechanism proposed in Q with the possibility of non-trivial flavour 
and momentum correlations, with initial- and final-state showers for all interactions, and with several 
options for colour correlations between initiator and remnant partons. Many of these improvements rely 
on the development of junction fragmentation in |4]. 

This is not the end of the line. Rather we see that many issues remain to understand better, such 
as colour correlations between partons in interactions and beam remnants, whereas others have not yet 
been studied seriously, such as the extent to which two interacting partons stem from the same initiator. 
Theoretical advances alone cannot solve all problems; guidance will have to come from experimental 
information. The increased interest in such studies bodes well for the future. 
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